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Abstract. We study thermal states of transiently accreting neutron stars (with mean accretion rates M ~ 10~^^ — 
10^" Mq yr~^) determined by the deep crustal heating of accreted matter sinking into stellar interiors. We 
formalize a direct correspondence of this problem to the problem of cooling neutron stars. Using a simple toy model 
we analyze the most important factors which affect the thermal states of accreting stars: a strong superfluidity in 
the cores of low-mass stars and a fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark 
phases of dense matter) in the cores of high-mass stars. We briefly compare the results with the observations of 
soft X-ray transients in quiescence. If the upper limit on the quiescent thermal luminosity of the neutron star 
in SAX J1808. 4-3658 (Campana et al. I2002II is associated with the deep crustal heating, it favors the model of 
nucleon neutron-star cores with switched-on direct Urea process. 
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1. Introduction 

Great progress in observations of soft X-ray transients 
(SXRTs) in quiescence has attracted attention to these 
objects. We consider the SXRTs containing neutron stars 
(NSs) in binary systems with low-mass companions (low- 
mass X-ray binaries); see Chen et al. p997f) for a re- 
view. They undergo the periods of outburst activity (days- 
months, sometimes years) superimposed with the periods 
of quiescence (months-decades). This transient activity is 
regulated most probably by the regime of accretion from 
the disks around the NSs. An active period is associated 
with a switched-on accretion; the accretion energy released 
at the NS surface is high enough for a transient to look 



like a bright X-ray source {Lx 
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The accretion is switched off or strongly suppressed dur- 
ing quiescence intervals when the NS luminosity drops by 
several orders of magnitude {Lx ^ 10^^ erg s~^). 

The nature of the quiescent emission is still a subject 
of debates. The hypothesis that this emission is produced 
by the thermal flux emergent from the NS interior has 
been rejected initially due to two reasons. First, the ra- 
diation spectra fitted by the blackbody model have given 
unreasonably small NS radii. Second, NSs in SXRTs have 
been expected to be internally cold; their quiescent emis- 
sion should have been much lower than the observed one. 
These arguments were questioned by Brown et al. H1998p . 
They suggested that the NSs can be warmed up to the 
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required level by the deep crustal heating associated with 
nuclear transformations in accreted matter sinking in the 
NS interiors (Haensel & Zdunik [1990(1 . while the radia- 
tion spectra can be fitted with the NS hydrogen atmo- 
sphere models (with realistic NS radii) . It turned out that 
the emergent radiation flux may depend on the NS in- 
ternal structure which opens an attractive possibility to 
explore the internal structure and the equation of state 
of dense matter by comparing the observations of SXRTs 
in quiescence with theoretical models (e.g., Ushomirsky & 
Rutledge HUUTl Colpi et al. I^UHTl Rutledge et al. li^OOi^al 
Brown et al. I2UU2I and references therein) . 

In this article, we discuss this possibility in more de- 
tail by making use of the close relationship between the 
theory of thermal states of transiently accreting NSs and 
the theory of cooling isolated NSs. Since the equation of 
state of NS matter is still poorly known we will make a 
general analysis of the problem with a simple toy model of 
NS thermal structure (described by Yakovlev & Haensel 
12003(1 . It will enable us to confront the theory with the ob- 
servations without performing complicated calculations. 

2. Numerical simulations 

Let us study a thermal state of a transiently accreting NS. 
Since the quiescence intervals are much shorter than typi- 
cal time scales of NS thermal relaxation 10'' yrs, Colpi 
et al. 12001(1 we neglect short-term variability in the NS 
crust; it can be associated with variable residual accre- 
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tion in quiescence, thermal relaxation of transient deep 
crustal heating, etc.; see, e.g., Ushomirsky & Rutledge 
1)2001(1 and Brown et al. H2002|l . Instead, we focus on a 
(quasi)stationary steady state of the NS determined by 
the mass accretion rate M = (M) (from 10~^^ to 10"^ 
Mq yr~^) averaged over time intervals comparable with 
the thermal relaxation time scales. The accretion rates 
of study are too low to noticeably increase NS mass, M, 
during long periods of the NS evolution. An accreted ma- 
terial sinks gradually into the deep layers of the NS crust 
(into the density range from - 10^° g cm~3 to - 10^^ g 
cm~^), where it undergoes various transformations (beta 
captures, neutron emissions and absorptions, pycnonu- 
clear reactions) accompanied by a substantial energy re- 
lease. This produces the deep crustal heating whose power 
is determined by M. Complete sequence of nuclear trans- 
formations and associated energy release was calculated 
by Haensel & Zdunik l)1990|l (see Bisnovatvi-Koean l^OOll 
for references to some earlier work). The total energy re- 
lease is about 1.45 MeV per accreting nucleon. The main 
energy release takes place at densities from about 10^^ to 
10^'^ g cm~'^, about 1 km under the surface. The heating 
power is estimated as 

Ldh = 1.45 MeVM/mN w 8.74 x lO^^ M^^ erg s"\ (1) 

where Mio = M/(10"^° Mg yr"^) and ton is the nucleon 
mass. Independently of the initial thermal state (cold or 
hot), a NS will eventually reach a stationary thermal state 
supported by this heating. 

We stress that Ldh is much smaller than the mean 
energy release rate at the NS surface (where kinetic ac- 
cretion energy transforms into heat) and in the outermost 
surface layers (at densities p ^ 10^° g cm~'^, where the ac- 
creted hydrogen or helium burns into heavier elements). 
We expect that the energy released at p ^ 10^° g cm~'^ 
is radiated away by the surface photon emission and does 
not heat the NS interior. This assumption is especially 
true for transiently accreting NSs, where the surface heat 
is rapidly emitted from the surface at the beginning of 
quiescent stages. 

The deep-heating theory of Haensel & Zdunik ( 199(|j) 
assumes that the burning in the very surface layers leads to 
the production of iron-like elements at p '--^ 10^° g cm~^. 
According to the recent studies by Schatz et al. (200l"| 
thermonuclear burning of hydrogen in the surface layers 
of accreting NSs may end at much heavier elements such 
as Sn, Sb and Te. The associated deep heating has been 
most recently analyzed by Haensel & Zdunik 1)2003(1 with 
the result that the heavy elements in the surface layers do 
not significantly affect the composition of accreted matter 
and the energy release in the deeper layers of the inner 
crust (where the pycnonuclear reactions proceed); the to- 
tal energy release remains about (1.1-1.5) MeV per nu- 
cleon. 

To study the thermal states of accreting NSs we em- 
ploy a simple toy model of coohng isolated NSs described 
by Yakovlev & Haensel 1(2003)1 . We have tested the results 



Table 1. Main processes of slow neutrino emission in nu- 
cleon matter: Murca and bremsstrahlung (brems) 



Process 






Qs, erg cm ^ s ^ 


Murca 


nN ^ 


pNei/ pNe nNu 


10^° - 3 X 10^^ 


Brems. 


NN ^ 




10" _ 1020 



with those obtained with the "exact" codes which calcu- 
late the thermal states of cooling NSs and accreting NSs; 
the agreement has turned out to be satisfactory. Taking 
into account many theoretical and observational uncer- 
tainties, the toy model seems sufficient at the present stage 
of the investigation. 

The toy model assumes that the NS core is divided 
into three zones: the outer zone, p < ps] the transition 
zone, Ps p < Pf, and the inner zone, p > pf. If the NS 
central density Pc < Ps, two last zones are absent. 

In the outer zone, the neutrino emission is supposed 
to be slow, while in the inner zone it is fast. The neutrino 
emissivity (erg s^^ cm^^) is assumed to be: 

"(P < Ps) = QsT|, Q'r\p>Pf)^QfTl (2) 

Here, Tg is the internal stellar temperature T expressed 
in 10^ K, while Qs and Qi are constants. For simplicity, 
the toy model uses the linear interpolation in p between 
(5t'°^ and Qjf^' in the transition zone. The approximation 
is crude but sufficient for exploring the main unknown 
features of the transition zone: its position and thickness. 

This generic description of covers many physical 
models of nucleon and exotic supranuclear matter with 
different leading neutrino processes listed in Tables 1 and 
2 (from Yakovlev & Haensel 2003 with the kind permission 
of the authors). In these tables, N is a nucleon (neutron or 
proton, n or p); e is an electron; v and 9 are neutrino and 
antineutrino; q is a quasinucleon (mixed n and p states); 
while u and d are quarks. 

In particular, Qs describes the modified Urea (Murca) 
process in a nonsuperfluid nucleon matter, or weaker NN- 
bremsstrahlung (e.g., nn-bremsstrahlung if Murca is sup- 
pressed by a strong proton superfluidity as considered 
by Kaminker et al. 12002(1 . The factor Qf describes the 
processes of fast neutrino emission: a very powerful di- 
rect Urea (Durca) process in nucleon matter (Lattimer 
et al. I1991() or somewhat weaker similar processes in ex- 
otic phases of matter (pion condensed, kaon condensed, 
or quark matter) as reviewed, e.g., by Pethick 1(1992(1 . If 
hyperons appear in the NS core in addition to nucleons, 
the neutrino emissivity is expected to be qualitatively the 
same as in the nucleon core. The bottom line of Table 2 
corresponds to a nonsuperfluid quark matter in NS cores. 

The transition zone models a switch-on of the fast neu- 
trino emission with growing p. In a nonsuperfluid matter, 
the lower density ps is a threshold density of the fast emis- 
sion; this threshold is usually sharp, i.e., Ps ~ pi- The su- 
perfluidity suppresses the neutrino emission. Even if the 
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Table 2. Leading processes of fast neutrino emission in 
nucleon matter and three models of exotic matter 



Model 


Process 






Qt, erg cm~^ s~^ 


Nucleon matter 


n 




pc - 


-+ nu 




Pion condensate 


q ^ 


qeu 


qe - 




10^3 _ iq26 


Kaon condensate 




qev 


qe - 


■i qv 


1023 _ IQ24. 


Quark matter 


d ^ 


uev 


ue - 


^ Av 





superfluidity strength is high at the formal (nonsuperfluid) 
threshold, it decreases with growing p in all realistic mi- 
crophysical models of superfluid matter. In this case, ps 
is the density where the superfluid suppression of the fast 
neutrino emission becomes weaker, while pf is the den- 
sity where this suppression almost vanishes. Thus the su- 
perfluidity can broaden a switch-on of the fast neutrino 
emission and produce a sufhciently wide transition zone. 
These effects are explained, e.g., by Yakovlev et al. (|20Q2|, 
taken nucleon Durca process and strong proton superflu- 
idity as an example. In addition to reducing the ordinary 
neutrino emission processes, the superfluidity initiates a 
specific neutrino emission associated with Cooper pairing 
of baryons (Flowers et al. il976jl . We neglect this emission 
in the present calculations but outline its consequences in 
Sect. 4. 

The toy model solves the equation of thermal balance 
of a cooling NS in the approximation of isothermal in- 
terior (see below) under a number of simplified assump- 
tions on NS internal structure (Yakovlev & IIaensel .2003 ) . 
The density profile in the star is approximate, p(r) = 
Pc (1 -rV-R^), so that the NS mass is M = 87tR^Pc/15, R 
being the stellar radius. The stiffness of various equations 
of state can be mimicked by choosing different M — R rela- 
tions. For simplicity, following Yakovlev & Haensel J2003'), 
we set R= 12 km and vary pc from 7 x 10^'* to 1.4 x 10^^ 
g cm~3, varying thus M from 1.02 Mq to 2.04 Mq. More 
realistic M — R relations will not change our principal 
conclusions. Various cooling regimes are regulated by four 
parameters: Qt, Qs, Ps, and pf, and are studied (Sect. 4) 
without complicated computation. 

We have updated the toy model by incorporating the 
deep crustal heating as described below. 



3. Relation to cooling of isolated stars 

A thermal state of an accreting NS is similar to a state of 
a cooling isolated NS with isothermal interior. 

Both states can be found by solving the equation of 
thermal balance of a NS with isothermal interior (e.g.. 
Glen & Sutherland [Mil: 



m 

dt 
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Fig. 1. Simultaneous calculation of cooling and heating 
curves for low-mass (thick lines) and high-mass (thin lines) 
NS models (see the text). Displayed are: the neutrino, pho- 
ton and total luminosities, L^, and L^^. (left vertical 
scale), of a cooling NS versus its age, and the accretion 
rate M (right vertical scale) of an accreting NS star with 
-^dh = -^tof Dotted horizontal line: Lf^^ for M = 10-^° 
Mq yr~^. Filled circles and dotted vertical arrows refer 
to those low-mass and high-mass cooling NSs which are 
equivalent to accreting NSs with the indicated M; open 
circles: for these cooling NSs. 



where Tg is the effective surface temperature, Ti{t) = 
r(r, t) e* is the redshifted internal temperature which is 
constant throughout the isothermal interior {p ^ Ph ^ 
10^" g cm^'^) with account for the effects of General 
Relativity; T(r, t) is the local internal temperature of mat- 
ter, and <I>(r) is the metric function. The relation between 
Ti and Tg is known from the solution of the thermal con- 
duction problem in the outermost heat-blanketing stellar 
envelope {p < ph)- Furthermore, C is the total heat ca- 
pacity of the star, = AnaT^ R^ (1 — r^/R) is the pho- 
ton surface luminosity as detected by a distant observer 
(rg = 2GM /c^ being the gravitational radius), and 



,dh 



47r 



2$ 



dr Qu,dh e 
^1 - 2Gm/{c'^r)' 



(4) 



(3) 



are the redshifted neutrino luminosity and the redshifted 
power of the deep crustal heating; Q^, is the neutrino 
emissivity; Qdh is the energy release rate due to the deep 
crustal heating, and m(r) is the gravitational mass con- 
tained in a sphere with radial coordinate r. In addi- 
tion, we introduce the effective surface temperature de- 
tected by a distant observer, = Tg \/l — r^/R, and 
the "apparent" NS radius, = i?/^l - r^/R. Then 

As mentioned above, we replace a transient accretion 
by a steady-state accretion with the mean constant ac- 
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Fig. 2. Heating curves of accreting NSs compared with 
observations of several SXRTs in quiescence. Thick sohd 
hne: deep crustal heating power, Ldh- Long dashes and 
thin sohd hnes: two models of heating curves of low-mass 
NSs (Qs = 3 X 10^^ and lO^^). Dotted, short-dashed, 
and dot-and-dashed lines: three models of high-mass NSs 
(Qf ^ 1023, 1025, and 1027). 

cretion rate Af — (M). We will study thus the steady 
states of NSs accreting at constant rates. A steady-state 
accretion in General Relativity is characterized by a con- 
stant mass accretion rate M which determines constant 
number of accreting baryons passing through a sphere of 
any radial coordinate r per unit time for a distant ob- 
server (e.g., Thorne & Zvtkow fl977|l . Gravitational red- 
shift in a thin NS crust may be regarded as constant. In 
this approximation, sufficient for our problem, we have 
constant e* = y^l — rg/R throughout the crust, and 
Lfi^ = idh rg/R, where Ldh is given by Eq. (jlj 

(neglecting the effects of General Relativity). In princi- 
ple, the gravitational redshift has to be included into L^. 
However, since M in SXRTs is determined with large un- 
certainties, it is premature to include such effects, and we 
set = Ldh in subsequent calculations. 

Following the cooling of an isolated NS we calculate 
a cooling curve, the dependence of the surface tempera- 
ture or the photon surface luminosity on NS age, T^{t) 
or L^{t). In about t ~ 100 yrs (when Tg°° drops below 
3x 10^ K) the NS interior becomes isothermal. Subsequent 
cooling can be considered as a sequence of quasistation- 
ary states. It is described by Eq. Q with = 0. The 
examples are presented in Fig. ^ described below. 

Following the evolution of an accreting NS from Eq. 
we would find that the NS, starting from an arbitrary 
initial thermal state, eventually reaches a stationary state 
supported by the deep crustal heating (whose power, L^, 
is given by Eq. QJ). The steady state is reached when the 
total NS luminosity, = L^f -I- L°°, is balanced by L^. 



transiently accreting neutron stars 

This steady state is described by Eq. © in the stationary 
limit: 

L^,(M) = L-(r,) + L-(T,), (5) 

where LJ^ is known once M is specified. Since we are 
interested in the steady state, it is sufficient to solve Eq. 

. The solutions give us a heating curve, the dependence 
of the surface temperature or the photon luminosity on the 
mean accretion rate, T^{M) or L^{M). For the values 
of M ^ 10~^ M0 yr^^, we obtain the thermal states with 
Tg ^ 3 X 10^ K. The examples are presented in Fig. |21 
discussed below. 

Comparing the equations of NS cooling and accretion 
heating we see that any steady state of a NS accreting at 
a rate M is equivalent to a state of a cooling NS in such a 
moment of time t when L^^(t) = L'^{M). In other words, 
calculating a coohng curve, we find L^^{t), L^{t) and 
equivalent M{t) for any moment of time. Considering t as 
a parameter, we immediately get a heating curve L'^{M). 

This procedure is illustrated in Fig. ^ It shows a si- 
multaneous solution of the cooling and heating problems 
with the toy model (Sect. 2) which calculates the cool- 
ing of isolated low-mass and high-mass NSs. The masses 
are 1.16 and 2.04 M© (pe = 8 x lO" and 1.4 x lO^^ g 
cm^3) xhc toy-model parameters are: Qs = erg 
cm~3 g-i and Qi — 102^ erg cm"^ (appropriate to 
the Murca and Durca processes in nonsuperfluid matter) , 
Ps — 8 X 10^^ g cm^3, and pi = 10^^ g cm"^. The figure 
shows the decrease of L^, L^, and (left vertical axis) 
with time for cooling stars. The right vertical axis presents 
M for accreting NSs with Ltot = -^dh- For instance, taking 
M — 10^^" M© yr~^ (horizontal dotted line) and locating 
the points, where Ltot = ^dh (filled circles and vertical 
dotted arrows), we find that the low-mass NS accreting 
with the given M is equivalent to the cooling NS of age 
< « 2 X 10** yrs, while the high-mass accreting NS is equiv- 
alent to the cooling NS of age t « 170 yrs. The thermal 
states of low-mass and high-mass accreting NSs are seen 
to be drastically different (the appropriate photon lumi- 
nosities are marked by open circles). 

It is clear that any steady state of an accreting NS is 
independent of the heat capacity of matter and the ther- 
mal conductivity in the isothermal NS interior, although 
the relaxation to this state does depend on these quanti- 
ties. Accordingly, superfluidity in NS interiors affects the 
steady state only by regulating the neutrino emission pro- 
cesses. 

Although a general analysis of Eq. ^ is complicated, 
we can mention two limiting regimes. 

(i) The photon emission regime is realized in cold 
enough NSs, where L"^ <C L^, i.e., the heat released in 
the deep crust is carried away by the thermal conductiv- 
ity to the surface and emitted in the form of the surface 
photon radiation. This regime is equivalent to the pho- 
ton emission stage of cooling NSs. In this regime, Eq. ||SJ) 
reduces to: L^^{AI) = L^{Tb), which immediately gives 

an estimate T^g = TJW^ K - 2Mj^/'* (for ~ 10 km). 
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Accordingly, the surface temperature is determined by the 
accretion rate, and it is independent of the internal struc- 
ture of the NS. The internal temperature can be found 
then using the Tg — relation. For typical surface gravities 
gs ^ 2 X 10^^ cm s^^ from the relation of Gudmundsson 
et al. ifTUS^ we have T,q = Ti/lO^ K - O.lTjj. Thus, 

T,9 ~ 0.4v/m^. 

(ii) The neutrino emission regime is realized in warmer 
NSs, where ^ L^, i.e., the heat released in the deep 
crust is spread by the thermal diffusion over the star 
and carried away by the neutrino emission. This regime 
is equivalent to the neutrino stage in cooling NSs. The 
thermal balance equation then reads: L|^(M) = LJf (Tj), 
which gives T^. The surface temperature can be found then 
using the Tg — Tj relation. In this regime, the surface tem- 
perature Ts does depend on the internal structure of the 
star, on the equation of state, the presence or absence of 
enhanced neutrino emission mechanism, and on the na- 
ture of this emission. For the assumptions implied in our 
toy model, we have two distinctly different cases. 

(ii.a) In low-mass NSs {pc < Ps) we have slow neutrino 
emission, Qf""^ = Q^Tf^. Introducing Q,2i = Qs/ClO^i 
erg s~^ cm~^), we have Ly ^ h x lO^^Tf^Qs2i erg s~^. 
Then we obtain T^g ~ 0.18 (Mio/Qs2i)^/^, and T^q - 
1.4 (Mio/Qs2i)^^^^; the surface temperature weakly de- 
pends on the accretion rate and the neutrino emission 
level Qs (as a result of the strong temperature depen- 
dence of L^). The comparison with case (i) shows that 
the transition from the photon emission regime to the neu- 

1/3 

trino emission regime takes place at Mio ~ 0.13/Qj,2i ^'^d 

Ts6 1.2/(3s2i^- 

(ii.b) In high-mass NSs {pc > pi) the neutrino emission 
is mainly produced by the fast mechanism, Q^^^^ = Qi^ig, 
in the inner zone (p > pt) of the NS core. Introducing 
Qf27 = Qt/(10^^ erg s~^ cm~^), we have L^, 5 x 
lO^^Tfg Qf27 erg s-\ ^ 0.01 {Mw/QmY^^ and T.g - 
0.3 (Mio/Qf27)^^^^; Ts is again a weak function of Al and 
Qf. A comparison with case (i) shows that the transition 
between the photon and neutrino emission regimes occurs 

at Mio - 10-V\/Qf27 and T^e - O.l/gJ/r - We have a 
much colder NS than in case (ii.a), and the neutrino emis- 
sion regime persists to much lower M. 

4. Results and discussion 

Our numerical calculations with the toy model confirm 
this qualitative analysis. The results are presented in Fig. 
|5] which shows the surface photon luminosity versus accre- 
tion rate for low-mass and high-mass NSs. 

The thick solid curve presents the deep heating power, 
L^, which is the upper limit of L"^ for any accreting 
source. 

Moving from top to bottom, the next two lines refer 
to low-mass NSs with two types of slow neutrino emission 
appropriate either to neutron-neutron bremsstrahlung in 
the NS cores with a strong proton superfluidity (Qs — 
3 X 10^^) or to the Murca process in nonsuperfluid cores 



(Qs — 10^^). The strong proton superfluidity damps the 
Murca process and enables us to obtain hotter NSs, just as 
in the theory of cooling NSs (e.g., Kaminker et al. I2002|l . 

The three next lines refer to high-mass NSs with three 
types of fast neutrino emission appropriate to (from top 
to bottom) kaon-condensed matter (Qf = lO^'^), pion- 
condensed matter (Qf = 10^^), or nucleon matter with 
open Durca process (Qf = 10^^). NSs with hyperons cores 
are expected to cool at about the same rate as NSs with 
nucleon cores. 

The heating curves of low-mass NSs provide the up- 
per limit of , while the curves of high-mass stars give 
the lower limit of , for any particular equation of state 
(EOS) of NS interiors (for any set of the four neutrino 
emission parameters, in our case). Accordingly, the upper 
limits are shaded below the curves, and the lower lim- 
its are shaded above the curves. As in the cooling theory 
(Yakovlev & Haensel IM75|l . the upper and lower heating 
curves are almost insensitive to the position and width of 
the transition zone if pi and Ps are located anywhere be- 
tween 8 X 10^^ and 1.2 x 10^*^ g cm^'^. The upper curve 
is determined by the parameter Qs, while the lower one is 
determined by Qf. 

Varying the NS mass from the lowest values to the 
highest we obtain a family of heating curves which fill in 
the space in Fig. [3 between the upper and lower curve for 
a given EOS of dense matter. The group of NSs whose 
heating curves lie essentially between the upper and lower 
curves will be called medium-mass stars. Their central 
density falls within the transition layer between the slow 
and fast neutrino emission zones (ps ^ P Pf, Sect. 2), 
and their mass range is sensitive to the position and width 
of this layer. The same situation occurs in the theory of 
cooling NSs (e.g., Yakovlev & Haensel l2003|l . Thus, for a 
given EOS of dense matter we obtain its own upper and 
lower heating curves, and intermediate heating curves of 
medium-mass stars. 

In analogy with the cooling theory, the existence of a 
representative class of medium-mass stars (a smooth tran- 
sition from the upper to the lower heating curves with 
increasing M) depends on the relative width of the tran- 
sition zone, Ap/ps = (pf — ps)/ps, and on the contrast be- 
tween the fast and slow neutrino emissivities, Qt/Qs- To 
ensure the existence of this NS class for a sharp emissivity 
contrast, Qt/Qs 3> 10"^, we need a rather wide transition 
zone, Ap/ps ^ 0.1. On the other hand, this class will be 
available even for a negligibly narrow zone if the emissivity 
contrast is lower, Qf/Qs ^ lO'^. 

These results can be confronted with the observations 
of SXRTs containing NSs. Although our theoretical toy 
model is oversimplified, in Fig. [5] we present an example 
of such an analysis for five SXRTs: Aql X-1, Cen X-4, 
4U 1608-522, KS 1731-26, and SAX J1808.4-3658. The 
data are rather uncertain. Thus we plot the observational 
points as thick dots. 

The parameters of the selected sources are collected in 
Table 01 The mean mass accretion rate M is evaluated as 
M — AM/ At, where AM is a total mass accreted over 
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Table 3. Parameters of NSs in SXRTs 



Source 


M, M© yr-i 


i^, erg s-i 


rr, eV 


R°°, km 


Reference 


Aql X-1 


1.0 X 10"^° 


5.3 X 10^^ 


113 


15.9 


Rutledee et al. (l2002bll : 


Table 6, fit 10 


Cen X-4 


1.4 X 10"^° 


3.1 X 10^^ 


75*) 


12.9 


Rutledge et al. i^M)T\ : 


Table 4 


4U 1608-522 


4.2 X 10""' 


4.1 X 10^^ 


170*) 


9.4 


Rutledge et al. (IB^ : 


Table 2 


KS 1731-260 


5.1 X 10"'' 


4.3 X 10^^ 


110 


15 


Wijnands et al. I^FM^ : 


Table 1, fit 2 


SAX J1808.4-3658 


5.0 X 10"^^ 


1.0 X 10^^ 






Campana et al. (I2002II : 


Sect. 2.2 



nonredshifted 



a representative period of time At. Both, AM and At, 
should include active and quiescent periods, although AM 
is mainly accumulated in the outburst states. In principle, 
we need M averaged over thermal relaxation time scales, 
~ 10'* yrs, while the observations provide us with sparse 
data over periods not longer than several decades. For Aql 
X-1, the mean M has been estimated by Rutledge et al. 
pinn|l (their Sect. 5, an estimate from the RXTE/ASM 
hght curve history). For Cen X-4 and 4U 1608-522, we 
obtain M from Table 9 in Chen et al. H1997|l . Following 
these authors, we take At = 8.67 yrs for 4U 1608-522, 
which is a frequently bursting source (AM is estimated 
for 6 outbursts in the period from 1970 to 1979). We take 
At = 33.16 yrs for Cen X-4 (with AM given for the only 
two outbursts in 1969 and 1979), adding the period from 
1979 till now when no outbursts were observed. Since the 
active states may be very rare for this source, it is safer to 
consider the obtained M as an upper limit. KS 1731-260 
recently (about 1.5 years ago) returned into quiescence 
after having actively accreted for « 11.5 years. For this 
source, we take At=11.5-f 1.5=13 yrs. We estimate the 
mass AM accreted during the long outburst state from 
the value of the mean outburst flux given by Rutledge 
et al. l|2002al Sect. 3.1). The estimation is made in the 
same manner as in Chen et al. 1)19971 Sect. 5.1.4). Since 
the recurrence time is unknown, our value of M is most 
probably an upper limit. Finally, M for SAX J1808.4- 
3658 was estimated by Bildstein & Chakrabarty H2001|l 
and Campana et al. H2002p . 

The values of in Fig. |3 are meant to refer to the 
quiescent thermal luminosity from the NS surfaces. For 
the first four sources, these values are obtained from the 
values of and R°° given in Table El The values of 
and R°° were evaluated by the authors cited in Table|21by 
fitting the observed spectra with the hydrogen atmosphere 
models. Note that the values of the surface temperature 
for Cen X-4 and 4U 1608-522 are nonredshifted. We have 
redshifted them assuming M = 1.4 M© and i? = 12 km. 
The spectrum of SAX J1808. 4-3658 in quiescence is well 
fitted by a power law, i.e., no surface thermal emission 
has been detected. The estimates of the upper limit of 
are model dependent and range from about 10'^° erg 
s-i to 2.5 X 10^1 erg s"* (Campana et al. lMT^ . We take 



L!^ = 10^* erg with the notice that the actual surface 
luminosity may be much lower. 

If the interpretation of quiescent emission as the ther- 
mal emission from the NS surfaces is correct, then all five 
NSs are heated to the neutrino emission stage (L^ > 
L^)- Since is reliably determined by the theory 

(Haensel & Zdunik ITMil IW)^ for a known M, and 
is measured, one can immediately estimate the neutrino 
luminosity of any source from the thermal balance, Eq. 
Q: — — L"^. In all our cases is comparable 
with (Fig.EJ. 

As seen from Fig. El we can treat NSs in 4U 1608-52 
and Aql X-1 as low-mass NSs with very weak neutrino 
emission from their cores (suppressed by strong nucleon 
superfluidity). The NSs in Cen X-4 and SAX J1808.4-3658 
seem to require the enhanced neutrino emission and are 
thus more massive. The status of the NS in KS 1731-26 is 
less certain because of poorly determined M (see above). 
If the real value of M is close to that in TableElit may also 
require some enhanced neutrino emission. Similar conclu- 
sions have been made by several authors (particularly, by 
Usho mirsky fc Rut l edge IMITl Colpi et al . 1201711 Rutledge 
et al. nnnU mmal I2002bl Brown et al. Wijnands 
et al. 12002(1 with respect to some of these sources or se- 
lected groups. Colpi et al. (pOOl) presented also the heat- 
ing curves for specific models of low-mass and high-mass 
NSs with superfluid nucleon cores and suggested that by 
tuning nucleon superfluidity and NS masses one can ex- 
plain the data. Using the toy-model, we can present a 
general analysis of the problem for different EOSs of NS 
interiors (assuming, of course, that all the sources have to 
be interpreted in terms of one EOS). In this way we can 
quantify the assumptions on enhanced neutrino emission 
in terms of pion condensed, kaon condensed, and Durca- 
allowed nucleon models of matter. 

Disregarding the SAX source for the moment, we can 
treat the NS in Cen X-4 either as a high-mass NS (with a 
kaon-condensed or quark core) or as a medium-mass NS 
(with a pion-condensed, quark, or Durca-allowed nucleon 
core) ; thus we cannot determine the nature of superdense 
matter. If the data on SAX J1808. 4-3658 are really rel- 
evant for our analysis we have the only choice to treat 
the NS as a high-mass NS with the nucleon core (and the 
NS in Cen X-4 as the medium-mass NS with the nucleon 
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core). This would mean that NS cores do not contain ex- 
otic phases of matter. 

Our toy model is too flexible and does not allow us to 
fix the position of the transition layer in the stellar cores 
where the slow neutrino emission transforms into the fast 
one. Adopting a specific EOS of NS interiors (with this 
position determined by microphysics input) we would be 
able to construct the sequences of heating curves for the 
stars with different M, and attribute certain values of M 
to any source ("weigh" NSs in SXRTs, as proposed by 
Colpi et al. 120011 just as in the case of cooling isolated 
NSs considered, e.g., by Kaminker et al. l2002|l . We intend 
to do this in our future publications, using an exact cooling 
code and taking into account some effects neglected in our 
simplified approach. 

In particular, we will account for the presence of light 
elements on the NS surfaces: they change the thermal con- 
ductivity of the NS heat-blanketing envelope and the re- 
lation between the surface and internal temperature of 
NSs. The effect is well known for cooling NSs (Potekhin 
et al. ll997|l and has been applied recently to accreting NSs 
(Brown et al. 12001 • 

We will also carefully treat the effects of baryon su- 
perfluidity in NS interiors (just as for cooling NSs, see, 
e.g., Kaminker et al. 12002(1 . Particularly, we will study the 
effects of Cooper-pairing neutrino emission of baryons ne- 
glected in the present analysis. Under certain conditions, 
this neutrino emission would violate our general assump- 
tions on the neutrino emissivity Qu {T, p) and complicate 
our study. For instance, according to our estimates, the 
■^P2 neutron superfluidity with the maximum values of the 
density dependent critical temperature Tcn{p) from ~ 10* 
K to ~ 2 X 10^ K in the nucleon NS core would produce 
a powerful Cooper pairing neutrino emission and strongly 
affect the thermal states of accreting NSs. However, the 
same effect would initiate a really fast cooling of not too 
massive isolated NSs in contradiction with the observa- 
tions of isolated cooling NSs (e.g., Kaminker et al. I2002S). 
Thus the presence of the indicated neutron superfluidity 
can be rejected on these grounds. 

Because of the similarity between the heating and cool- 
ing curves, the observations of cooling isolated NSs and ac- 
creting NSs in SXRTs can be analyzed together employing 
the same EOSs of NS interiors. This increases the statis- 
tics of the sources and the confidence of the results. The 
theory of cooling NSs has recently been confronted with 
observations by Yakovlev & Haensel (2003). Some cooling 
NSs (first of all, RX J0822-43 and PSR 1055-52) can be 
interpreted (Kaminker et al. '2002) as low-mass NSs with 
strong proton superfluidity in their cores. Other sources 
(first of all. Vela and Geminga) seem to require a fast neu- 
trino emission but the nature of this emission (a choice of 
fast-cooling model from Table 2) is uncertain, just as for 
SXRTs disregarding the data on SAX J1808.4-3658. In 
this context, the latter source is now the only one which 
indicates the absence of exotic phases of matter in NS 
cores. 



The assumption that the observed X-ray emission of 
SXRTs in quiescence emerges from the NS interior is still 
an attractive hypothesis. In any case the theory of deep 
crustal heating is solid and leaves no doubts that this heat- 
ing does occur in accreting NSs leading to observational 
consequences. 
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